\(\int \csc ^3(e+f x) (a+a \csc (e+f x))^m \, dx\) [30]

   Optimal result
   Rubi [A] (verified)
   Mathematica [A] (verified)
   Maple [F]
   Fricas [F]
   Sympy [F]
   Maxima [F]
   Giac [F]
   Mupad [F(-1)]

Optimal result

Integrand size = 21, antiderivative size = 156 \[ \int \csc ^3(e+f x) (a+a \csc (e+f x))^m \, dx=\frac {\cot (e+f x) (a+a \csc (e+f x))^m}{f \left (2+3 m+m^2\right )}-\frac {\cot (e+f x) (a+a \csc (e+f x))^{1+m}}{a f (2+m)}-\frac {2^{\frac {1}{2}+m} \left (1+m+m^2\right ) \cot (e+f x) (1+\csc (e+f x))^{-\frac {1}{2}-m} (a+a \csc (e+f x))^m \operatorname {Hypergeometric2F1}\left (\frac {1}{2},\frac {1}{2}-m,\frac {3}{2},\frac {1}{2} (1-\csc (e+f x))\right )}{f (1+m) (2+m)} \]

[Out]

cot(f*x+e)*(a+a*csc(f*x+e))^m/f/(m^2+3*m+2)-cot(f*x+e)*(a+a*csc(f*x+e))^(1+m)/a/f/(2+m)-2^(1/2+m)*(m^2+m+1)*co
t(f*x+e)*(1+csc(f*x+e))^(-1/2-m)*(a+a*csc(f*x+e))^m*hypergeom([1/2, 1/2-m],[3/2],1/2-1/2*csc(f*x+e))/f/(m^2+3*
m+2)

Rubi [A] (verified)

Time = 0.24 (sec) , antiderivative size = 156, normalized size of antiderivative = 1.00, number of steps used = 5, number of rules used = 5, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.238, Rules used = {3885, 4086, 3913, 3912, 71} \[ \int \csc ^3(e+f x) (a+a \csc (e+f x))^m \, dx=-\frac {2^{m+\frac {1}{2}} \left (m^2+m+1\right ) \cot (e+f x) (\csc (e+f x)+1)^{-m-\frac {1}{2}} (a \csc (e+f x)+a)^m \operatorname {Hypergeometric2F1}\left (\frac {1}{2},\frac {1}{2}-m,\frac {3}{2},\frac {1}{2} (1-\csc (e+f x))\right )}{f (m+1) (m+2)}+\frac {\cot (e+f x) (a \csc (e+f x)+a)^m}{f \left (m^2+3 m+2\right )}-\frac {\cot (e+f x) (a \csc (e+f x)+a)^{m+1}}{a f (m+2)} \]

[In]

Int[Csc[e + f*x]^3*(a + a*Csc[e + f*x])^m,x]

[Out]

(Cot[e + f*x]*(a + a*Csc[e + f*x])^m)/(f*(2 + 3*m + m^2)) - (Cot[e + f*x]*(a + a*Csc[e + f*x])^(1 + m))/(a*f*(
2 + m)) - (2^(1/2 + m)*(1 + m + m^2)*Cot[e + f*x]*(1 + Csc[e + f*x])^(-1/2 - m)*(a + a*Csc[e + f*x])^m*Hyperge
ometric2F1[1/2, 1/2 - m, 3/2, (1 - Csc[e + f*x])/2])/(f*(1 + m)*(2 + m))

Rule 71

Int[((a_) + (b_.)*(x_))^(m_)*((c_) + (d_.)*(x_))^(n_), x_Symbol] :> Simp[((a + b*x)^(m + 1)/(b*(m + 1)*(b/(b*c
 - a*d))^n))*Hypergeometric2F1[-n, m + 1, m + 2, (-d)*((a + b*x)/(b*c - a*d))], x] /; FreeQ[{a, b, c, d, m, n}
, x] && NeQ[b*c - a*d, 0] &&  !IntegerQ[m] &&  !IntegerQ[n] && GtQ[b/(b*c - a*d), 0] && (RationalQ[m] ||  !(Ra
tionalQ[n] && GtQ[-d/(b*c - a*d), 0]))

Rule 3885

Int[csc[(e_.) + (f_.)*(x_)]^3*(csc[(e_.) + (f_.)*(x_)]*(b_.) + (a_))^(m_), x_Symbol] :> Simp[(-Cot[e + f*x])*(
(a + b*Csc[e + f*x])^(m + 1)/(b*f*(m + 2))), x] + Dist[1/(b*(m + 2)), Int[Csc[e + f*x]*(a + b*Csc[e + f*x])^m*
(b*(m + 1) - a*Csc[e + f*x]), x], x] /; FreeQ[{a, b, e, f, m}, x] && EqQ[a^2 - b^2, 0] &&  !LtQ[m, -2^(-1)]

Rule 3912

Int[(csc[(e_.) + (f_.)*(x_)]*(d_.))^(n_.)*(csc[(e_.) + (f_.)*(x_)]*(b_.) + (a_))^(m_), x_Symbol] :> Dist[a^2*d
*(Cot[e + f*x]/(f*Sqrt[a + b*Csc[e + f*x]]*Sqrt[a - b*Csc[e + f*x]])), Subst[Int[(d*x)^(n - 1)*((a + b*x)^(m -
 1/2)/Sqrt[a - b*x]), x], x, Csc[e + f*x]], x] /; FreeQ[{a, b, d, e, f, m, n}, x] && EqQ[a^2 - b^2, 0] &&  !In
tegerQ[m] && GtQ[a, 0]

Rule 3913

Int[(csc[(e_.) + (f_.)*(x_)]*(d_.))^(n_.)*(csc[(e_.) + (f_.)*(x_)]*(b_.) + (a_))^(m_), x_Symbol] :> Dist[a^Int
Part[m]*((a + b*Csc[e + f*x])^FracPart[m]/(1 + (b/a)*Csc[e + f*x])^FracPart[m]), Int[(1 + (b/a)*Csc[e + f*x])^
m*(d*Csc[e + f*x])^n, x], x] /; FreeQ[{a, b, d, e, f, m, n}, x] && EqQ[a^2 - b^2, 0] &&  !IntegerQ[m] &&  !GtQ
[a, 0]

Rule 4086

Int[csc[(e_.) + (f_.)*(x_)]*(csc[(e_.) + (f_.)*(x_)]*(b_.) + (a_))^(m_)*(csc[(e_.) + (f_.)*(x_)]*(B_.) + (A_))
, x_Symbol] :> Simp[(-B)*Cot[e + f*x]*((a + b*Csc[e + f*x])^m/(f*(m + 1))), x] + Dist[(a*B*m + A*b*(m + 1))/(b
*(m + 1)), Int[Csc[e + f*x]*(a + b*Csc[e + f*x])^m, x], x] /; FreeQ[{a, b, A, B, e, f, m}, x] && NeQ[A*b - a*B
, 0] && EqQ[a^2 - b^2, 0] && NeQ[a*B*m + A*b*(m + 1), 0] &&  !LtQ[m, -2^(-1)]

Rubi steps \begin{align*} \text {integral}& = -\frac {\cot (e+f x) (a+a \csc (e+f x))^{1+m}}{a f (2+m)}+\frac {\int \csc (e+f x) (a (1+m)-a \csc (e+f x)) (a+a \csc (e+f x))^m \, dx}{a (2+m)} \\ & = \frac {\cot (e+f x) (a+a \csc (e+f x))^m}{f \left (2+3 m+m^2\right )}-\frac {\cot (e+f x) (a+a \csc (e+f x))^{1+m}}{a f (2+m)}+\frac {\left (1+m+m^2\right ) \int \csc (e+f x) (a+a \csc (e+f x))^m \, dx}{(1+m) (2+m)} \\ & = \frac {\cot (e+f x) (a+a \csc (e+f x))^m}{f \left (2+3 m+m^2\right )}-\frac {\cot (e+f x) (a+a \csc (e+f x))^{1+m}}{a f (2+m)}+\frac {\left (\left (1+m+m^2\right ) (1+\csc (e+f x))^{-m} (a+a \csc (e+f x))^m\right ) \int \csc (e+f x) (1+\csc (e+f x))^m \, dx}{(1+m) (2+m)} \\ & = \frac {\cot (e+f x) (a+a \csc (e+f x))^m}{f \left (2+3 m+m^2\right )}-\frac {\cot (e+f x) (a+a \csc (e+f x))^{1+m}}{a f (2+m)}+\frac {\left (\left (1+m+m^2\right ) \cot (e+f x) (1+\csc (e+f x))^{-\frac {1}{2}-m} (a+a \csc (e+f x))^m\right ) \text {Subst}\left (\int \frac {(1+x)^{-\frac {1}{2}+m}}{\sqrt {1-x}} \, dx,x,\csc (e+f x)\right )}{f (1+m) (2+m) \sqrt {1-\csc (e+f x)}} \\ & = \frac {\cot (e+f x) (a+a \csc (e+f x))^m}{f \left (2+3 m+m^2\right )}-\frac {\cot (e+f x) (a+a \csc (e+f x))^{1+m}}{a f (2+m)}-\frac {2^{\frac {1}{2}+m} \left (1+m+m^2\right ) \cot (e+f x) (1+\csc (e+f x))^{-\frac {1}{2}-m} (a+a \csc (e+f x))^m \operatorname {Hypergeometric2F1}\left (\frac {1}{2},\frac {1}{2}-m,\frac {3}{2},\frac {1}{2} (1-\csc (e+f x))\right )}{f (1+m) (2+m)} \\ \end{align*}

Mathematica [A] (verified)

Time = 5.54 (sec) , antiderivative size = 178, normalized size of antiderivative = 1.14 \[ \int \csc ^3(e+f x) (a+a \csc (e+f x))^m \, dx=-\frac {(a (1+\csc (e+f x)))^m \left ((-2+m) m \cot ^4\left (\frac {1}{2} (e+f x)\right ) \operatorname {Hypergeometric2F1}\left (-2-m,-2 m,-1-m,-\tan \left (\frac {1}{2} (e+f x)\right )\right )+(2+m) \left (m \operatorname {Hypergeometric2F1}\left (2-m,-2 m,3-m,-\tan \left (\frac {1}{2} (e+f x)\right )\right )+2 (-2+m) \cot ^2\left (\frac {1}{2} (e+f x)\right ) \operatorname {Hypergeometric2F1}\left (-2 m,-m,1-m,-\tan \left (\frac {1}{2} (e+f x)\right )\right )\right )\right ) \tan ^2\left (\frac {1}{2} (e+f x)\right ) \left (1+\tan \left (\frac {1}{2} (e+f x)\right )\right )^{-2 m}}{4 f (-2+m) m (2+m)} \]

[In]

Integrate[Csc[e + f*x]^3*(a + a*Csc[e + f*x])^m,x]

[Out]

-1/4*((a*(1 + Csc[e + f*x]))^m*((-2 + m)*m*Cot[(e + f*x)/2]^4*Hypergeometric2F1[-2 - m, -2*m, -1 - m, -Tan[(e
+ f*x)/2]] + (2 + m)*(m*Hypergeometric2F1[2 - m, -2*m, 3 - m, -Tan[(e + f*x)/2]] + 2*(-2 + m)*Cot[(e + f*x)/2]
^2*Hypergeometric2F1[-2*m, -m, 1 - m, -Tan[(e + f*x)/2]]))*Tan[(e + f*x)/2]^2)/(f*(-2 + m)*m*(2 + m)*(1 + Tan[
(e + f*x)/2])^(2*m))

Maple [F]

\[\int \csc \left (f x +e \right )^{3} \left (a +a \csc \left (f x +e \right )\right )^{m}d x\]

[In]

int(csc(f*x+e)^3*(a+a*csc(f*x+e))^m,x)

[Out]

int(csc(f*x+e)^3*(a+a*csc(f*x+e))^m,x)

Fricas [F]

\[ \int \csc ^3(e+f x) (a+a \csc (e+f x))^m \, dx=\int { {\left (a \csc \left (f x + e\right ) + a\right )}^{m} \csc \left (f x + e\right )^{3} \,d x } \]

[In]

integrate(csc(f*x+e)^3*(a+a*csc(f*x+e))^m,x, algorithm="fricas")

[Out]

integral((a*csc(f*x + e) + a)^m*csc(f*x + e)^3, x)

Sympy [F]

\[ \int \csc ^3(e+f x) (a+a \csc (e+f x))^m \, dx=\int \left (a \left (\csc {\left (e + f x \right )} + 1\right )\right )^{m} \csc ^{3}{\left (e + f x \right )}\, dx \]

[In]

integrate(csc(f*x+e)**3*(a+a*csc(f*x+e))**m,x)

[Out]

Integral((a*(csc(e + f*x) + 1))**m*csc(e + f*x)**3, x)

Maxima [F]

\[ \int \csc ^3(e+f x) (a+a \csc (e+f x))^m \, dx=\int { {\left (a \csc \left (f x + e\right ) + a\right )}^{m} \csc \left (f x + e\right )^{3} \,d x } \]

[In]

integrate(csc(f*x+e)^3*(a+a*csc(f*x+e))^m,x, algorithm="maxima")

[Out]

integrate((a*csc(f*x + e) + a)^m*csc(f*x + e)^3, x)

Giac [F]

\[ \int \csc ^3(e+f x) (a+a \csc (e+f x))^m \, dx=\int { {\left (a \csc \left (f x + e\right ) + a\right )}^{m} \csc \left (f x + e\right )^{3} \,d x } \]

[In]

integrate(csc(f*x+e)^3*(a+a*csc(f*x+e))^m,x, algorithm="giac")

[Out]

integrate((a*csc(f*x + e) + a)^m*csc(f*x + e)^3, x)

Mupad [F(-1)]

Timed out. \[ \int \csc ^3(e+f x) (a+a \csc (e+f x))^m \, dx=\int \frac {{\left (a+\frac {a}{\sin \left (e+f\,x\right )}\right )}^m}{{\sin \left (e+f\,x\right )}^3} \,d x \]

[In]

int((a + a/sin(e + f*x))^m/sin(e + f*x)^3,x)

[Out]

int((a + a/sin(e + f*x))^m/sin(e + f*x)^3, x)